"Phase Diagram" of the Spin Hall Effect 
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We obtain analytic formulas for the frequency-dependent spin-Hall conductivity of a two- 
dimensional electron gas (2DEG) in the presence of impurities, linear spin-orbit Rashba interaction, 
and external magnetic field perpendicular to the 2DEG. We show how different mechanisms (skew- 
scattering, side-jump, and spin precession) can be brought in or out of focus by changing controllable 
parameters such as frequency, magnetic field, and temperature. We find, in particular, that the d.c. 
spin Hall conductivity vanishes in the absence of a magnetic field, while a magnetic field restores 
the skew-scattering and side-jump contributions proportionally to the ratio of magnetic and Rashba 
fields. 
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The spin-Hall effect (SHE), i.e. the generation of a 
transverse spin current in response to a d.c. electric 
h>k]l— — — ~ has attracted much attention recently, par- 
ticularly after a series of experimentsZ'&S'iS culminating 
in the observation of the SHE at room temperature^. 
On the theoretical front, however, there remains consid- 
erable uncertainty as to the physical origin of the SHE, 
which appears to depend on an intricate interplay of three 
processes: (i) the skew-scattering (SS) due to spin-orbit 
interaction (SOI) between electrons and impurities, (ii) 
the side-jump (SJ) (due to the non-canonical character of 
the physical position and velocity variables in the pres- 
ence of SOI with impurities), and (iii) the spin precession 
caused by spin non-conserving terms in the band struc- 
ture, among which we include the linear Rashba SOI gen- 
erated by an external electric field. To these we may add 
the influence of an external magnetic field, which tends 
to lock the spins in a fixed direction, thus reducing the 
importance of spin precession. A first principles theory 
that include all of these effects on equal footing is very 
complicated. To our knowledge, the diagrammatic ap- 
proach by Tse and Das Sarma— comes closest to fulfill- 
ing the order, and yet it does not include magnetic field 
or frequency. However, these diagrammatic calculations 
are very difficult to follow in detail and do not lead to an 
intuitive understanding of the striking nonadditivc be- 
havior of impurity and band structure effects. 

Our goal in this paper is to present the "phase dia- 
gram" of the SHE, i.e. to clarify in which range of ex- 
perimentally controllable parameters one should expect 
the dominance of each mechanism mentioned above, and 
how the crossovers between different regimes occur. We 
do this for a 2DEG with Rashba SOI and a magnetic 
field perpendicular to the plane. The two parameters 
that are most easily controlled in an experiment are (i) 
the frequency u) of the electric field and (ii) the magni- 
tude of the perpendicular magnetic field, whose strength 
is characterized by the resonance frequency luq. Accord- 
ingly, we plot our "phase-diagram" in the lot — ujq/ a^kp 
plane, where r is the electron impurity scattering time 
and ankp is the magnitude of the effective magnetic field 



due to Rashba spin-orbit coupling for electrons at the 
Fermi wavevector kp. Throughout the paper we assume 
u), 1/t, a^kp, Wo Ep, i.e. all energy scales are smaller 
than the Fermi energy. 

Our qualitative conclusions are shown in Fig. ([!]). In 



Skew scattering/ 
side jump 

low frequency and 
strong magnetic field 



or 

o 

3 



~1 



low frequency and 
weak magnetic field 



a - ml (-2a p 3 \ 

KM \ M 2 1 / 




KM 



high frequency and 
weak magnetic field 



0=0 



-1 



COT 



FIG. 1: Different regimes of spin-Hall effect in the cut — 
uio/aRkF plane, o\d is Drude conductivity, p ss is skew scat- 
tering. Spin-orbit interactions and mobility are fixed. See 
discussion in the text. 

contrast to previous calculations 11 , we found that the d.c. 
limit (ur — > 0) of the spin Hall conductivity is zero in 
the presence of spin precession and in the absence of a 
magnetic field (u>q = 0). As the magnetic field increases, 
both SS and SJ contributions increase with the ratio of 
the magnetic to the Rashba field (uj /ak F ), thus restor- 
ing the values they would have had at zero frequency 
and zero magnetic field in the absence of spin preces- 
sion. In high-mobility samples the SS mechanism is the 
dominant mechanis m 12 i 13 overcoming the SJ contribu- 
tion, which has an opposite sign for attractive impurity 
potential. However, the SJ mechanism could well domi- 
nate in low-mobility samples. As discussed in Ref.— , the 
mobility can be controlled to some extent by changing 
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the temperature T, and this allows one in principle to 
tune between SS and SJ contributions. For this reason 
we use the label SS/SJ for the left side of diagram, where 
lot <C 1. In the a.c. regime (lot 3> 1) and for low mag- 
netic field impurities become irrelevant, leaving room for 
the intrinsic SHE with "universal" conductivity e/8ir.— 
Finally, in the a.c. regime and at high magnetic field the 
spin Hall conductivity declines to zero in different man- 
ners depending on whether lot or the magnetic field is 
kept constant as shown in Fig. |T]). 

In what follows we describe the main points of our 
new theoretical approach, which enables us to calculate 
the spin Hall conductivity in different regions of the pa- 
rameter space and to derive simple analytic formulas de- 
scribing the crossovers between different regimes. 

Our model is defined by the hamiltonian 
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Besides the kinetic energy and the usual electron- 
impurity potential V(r), we have two distinct spin-orbit 
couplings. The a\ coupling between the electrons and the 
impurities conserves the z-component of the spin and is 
responsible for SS and the SJ effects. The an coupling 
also known as Rashba coupling - creates a momentum- 
dependent magnetic field in the plane, which breaks the 
conservation of S z and causes spin precession. This term 
is responsible, under appropriate conditions, for the in- 
trinsic contribution to SHE. Finally, we have included a 
magnetic field perpendicular to the plane. The Zeeman 
splitting, wo, and the frequency lo of the a.c. electric field 
(not shown in H ) are the two control knobs in terms of 
which our "phase diagram" will be plotted. 

The skew-scattering effect is easily described in the 
Boltzmann equation formalism^, but it is difficult to 
treat in the diagrammatic approach^. On the other 
hand, spin precession effects are easily included in the di- 
agrammatic formalism, but are problematic in the Boltz- 
mann equation formalism (the distribution function be- 
comes a 2 x 2 matrix) . We get the best of two worlds by 
combining the two approaches in the following manner. 
First we notice that the skew-scattering collision term in 
the Boltzmann equation is formally equivalent to the im- 
position of a "spin-electric field" E y , which accelerates 
up-spin and down-spin electrons ("up" and "down" are 
defined with respect to a z-axis) in opposite directions 
along the y axis, perpendicular to the flow of the charge 
current (x). The problem is now "reduced" to calculating 
the z-spin current j y which flows along the y axis in re- 
sponse to the spin electric field E y in the same direction. 
This can be done with the help of the standard diagram- 
matic formalism, including both electron-impurity scat- 
tering and spin precession, but not the skew-scattering 
processes; for the skew scattering has already "done its 
job" by producing the spin electric field Ef. (in other 



words, we work to first order in a±). Thus we have 

f y \ ss = <j z yy E* (2) 

where o~ yy is the longitudinal spin conductivity calculated 
in the absence of skew-scattering (or side-jump) effects. 
On the other hand, E y has the well-known expression 

= P SS j X (3) 

where j x is the current density in the x direction and 
p ss = —m/ne 2 T ss is the skew-scattering resistivity, calcu- 
lated from the Boltzmann collision integral^ An explicit 
expression for the skew-scattering rate 1/t ss is given in 
Eq.(29) of Ref. 12 : notice that it is proportional to ct\ and 
positive. 16 - j x can be expressed as arjE x , where E x is the 
electric field in the x direction and ou is the Drude con- 
ductivity of the electron gas. Therefore Eq. [2] becomes: 

fyU = a z yy p ss a D E x , (4) 

from which we extract the first important result of this 
paper, namely the expression for the skew-scattering con- 
tribution to the spin Hall conductivity a SH , 

= <TyyP SS 0- D (5) 

The dependence of this formula on spin-precession rate, 
frequency of the a.c. field, and magnetic field will be 
obtained below. In particular we will show that crf/^ 
vanishes at low magnetic field (spin-precession regime), 
and it recovers the zero-precession value (o~ 2 D p ss h/ e) at 
high magnetic field. 

We will then consider separately the remaining con- 
tributions afP and o-p H to the spin Hall conductivity. 
They are most efficiently analyzed in terms of the Kubo 
formula for the spin-Hall conductivity. The SJ contribu- 
tion, similar to SS term vanishes at low magnetic field 
(spin- precession regime), and recovers its value in high 
magnetic fields. The remaining intrinsic contribution is 
easily calculated by standard diagrammatics (including 
vertex corrections), and leads to the well known e/8ir 
result in the appropriate a.c. regime. 

Skew-scattering - As discussed above, the central role 
in calculating aff 1 is played by the longitudinal spin- 
channel conductivity a y . The formal expression for the 
real part of a yy is 

(J yy(u) = -7TT (t>) 

m n lo 

where p y , S z are momentum and spin operators for a sin- 
gle electronic, the double bracket is the usual notation 
for the spin-current-spin current response function, n is 
the areal density of the electron gas. Our objective is 
to calculate this Kubo formula including elastic electron- 
impurity scattering in the Born approximation, spin pre- 
cession due to a Rashba spin-orbit interaction and a mag- 
netic field along the z-axis, but neglecting spin-orbit in- 
teractions with the impurities. The Feynman diagrams 
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FIG. 2: Diagrammatic representation of ladder approxima- 
tion to calculate the spin-channel conductivity. 



for <jy y {uj), to the desired level of accuracy, are shown in 
Fig. 2. The solid lines represent 2x2 Green's functions 
including Rashba spin-orbit coupling, magnetic field, and 
a diagonal disorder self-energy (i/2r)sgn(oj), where 1/r 
is the elastic scattering rate: 



G(p,w) 



lu - e p + h p S + ^rsgn(cj) 
(^-£ p +27SgnM) 2 -/i2/4 



(7) 



Here h p = (otRk y , — aRk x -, too) combines the Rashba in- 
plane field (x- and y-components) and the magnetic field 
along the z axis in a single effective field. The dashed 
lines are electron-impurity interactions (averaged over 
disorder), and there are spin-current operators p y S z at 
the vertices. The empty bubble, calculated by the stan- 
dard procedure with ctRkp, loq, and 1/r all much smaller 
than the Fermi energy is 
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(8) 



where ctd = ne 2 T/m* is the Drude conductivity and 
fl = \J (aflfcp) 2 + is the spin precession frequency. 
Eq. {S| suggests that the spin channel conductivity is 
similar to the ordinary Drude conductivity, only some- 
what reduced by precession about the in-plane Rashba 
field. The inclusion of vertex corrections changes the re- 
sult drastically. The sum of the ladder diagrams produces 
a "dressed" spin-current vertex A y of the form: 
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aRk 2 F [(aj T)S x - (1 - iuJT)Sy] 
a 2 R k 2 F r - 2VL 2 t{1 - iur) 



(9) 



Including A y in the calculation of the spin-channel con- 
ductivity we get 



2h<j D cos 2 6(1 + cos 2 S) 
e (1 + cos 2 S) 2 + Alj 2 t 2 



(10) 



where cos S = ^ is the cosine of the angle between the 
external magnetic field and the full effective magnetic 
field Q. While Eq. (fT0|) reduces to the empty-bubble 
result in the limit loo 3> ctRkF, its most striking feature 
is that it vanishes identically (i.e. at all frequencies) in 
the opposite limit loq <C ctRkF- In other words, the spin 
channel conductivity is zero in the absence of an external 
magnetic field, as long as an infinitesimal spin precession 
rate chrIsf is present! Accordingly, the skew scattering 
spin-Hall conductivity from Eq. ([5]) is^ 



SH 2hu 2 D p ss cos 2 (5(1 + cos 2 6) 
Tss = e (1 + cos 2 S) 2 + 4uj 2 t 2 



(11) 



which vanishes in the absence of an external magnetic 
field and recovers the value a 2 D p ss h/e for strong mag- 
netic fields. The vanishing of the spin-conductivity is a 
peculiar feature of the linear Rashba model, in which the 
spin current p y S z is proportional to the time derivative 

of S y : p y S z — —Sy/aR. The expectation value of a time 
derivative must vanish at zero frequency. This is the 
same reason which causes the vanishing of the Rashba 
spin-Hall conductivity in the d.c. limit. 

Side-jump - In a recent paper— we have shown how to 
identify the SJ contribution in the Kubo linear response 
formalism. We perturb the hamiltonian (fTJ) with a uni- 
form electric field of frequency uj in the x direction. Then 
after a series of manipulations, which have been described 
in Ref. — , we arrive at the following Kubo formula for 
the spin Hall conductivity: 
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The first term on the right hand side produces SS as well 
as possible "intrinsic" contributions which are discussed 
in next section. The second line of this equation, which 
explicitly shows a dependence on a\, is responsible for 
the SJ contribution 19 , denoted a^ n H . We now focus on 
this contribution. Using the Heisenberg equation of mo- 
tion for the momentum operator pi x = —\/ x V(fi) in zero 
order in ai (because SJ terms are already explicitly linear 
in ax), we rewrite the SJ contribution as 



ea\n 
2im*uj 



((Px;Px)) ~ 



mpyS z ;p y S z )) 

h 2 



(13) 



In the first term in the square brackets we apply the 
standard rule of integration, which allows us to replace 
Px/iu by —fix and rewrite it (its real part) as: 

ceien ■ —ct\ne A aine 

;{{Px;Px)) = 7T-1 ((Px,Px)) = Wf. ; ; , 2 ox , (14) 
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2m* 



2(1 +UJ 2 T 2 ) 



where the last equality follows from the well-known form 
of the current-current response function in a weakly dis- 
ordered system. The second term in the square brackets 
of Eq. (|13[) can be rewritten as follows: 
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/vS i>,,S = -^-((PySzl 4 (PySz) -PySz)) (15) 



Applying again the integration formula, the first term of 
Eq. El yields: 

— 2enu\ llr . - „ a ect\n cos 2 <5(1 + cos 2 5) 2 
—-T—((PyS z ;p y S z )) - 



2 (1 + cos 2 <5) 2 +4w 2 t 2 



(16) 

In the absence of spin precession (i.e., for lu » a^fci?, 
cos5 ~ 1) and for zero frequency the contributions (I14p . 
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which in the d.c. limit simplifies to 




Sfl, 



FIG. 3: Behavior of spin-Hall conductivity as a function of lot 
and LJo/ankF and carrier concentration ri2D — 2 x 10 12 cm -2 , 
fi =lm 2 /Vs, m = 0.067m e , ai =0.053nm 2 ,a B fc F =10meV. 
For calculation of skew-scattering, we assumed square poten- 
tial characterized by t/t ss = 0.002 like in Ref^ 
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(18) 



Summarizing, the side-jump conductivity in the d.c. 
limit (wt <C 1) grows from zero at uj — to the "full" 
value (aine) at high magnetic fields, while it tends to 
zero in high frequency limit. 

Rashba contribution - Evaluating the first line of 
Eq. (fT2)l in the absence of skew-scattering and side- 
jump, which we have already taken into account, leads to 
the well-known resul t 14 ' 15 a^ H (oj) = j— ^L^ ) which 
vanishes at to — and tends to the "ballistic" limit e/87r 
for lu ^> 1/t. We have found that the behavior of af> H (uj) 
as a function of external perpendicular magnetic field and 
frequency is given by: 



SH 



(I16p add up to the "usual" SJ conductivity a\ne of 
RefJ^. In the opposite limit of u>q <C ot^kp (cosS = 0) 
this term vanishes for the same reason as o~ s . The 
second term on the rhs of Eq. (|15p can be rewritten 

as ({S y ;p y (p ■ S))} (we use the equations of motion 

S z = ctRP ■ S and p y S z — —Sy/an). In the limit of 
zero frequency and zero magnetic field this reduces to 
{(pySy]p y Sy)) = ( H 2 /4) { (p x ; p x ) ) , which exactly cancels 
the contribution from (|14[) and causes the S J contribution 
to vanish. For general magnetic fields and frequencies the 
total afj* is finally given by the formula 2 ^ 
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ot\ne 



i 2 S 



2cos 2 ^(l + cos 2 S) 
1 + w 2 t 2 ^ (1 + cos 2 5) 2 + 4lu 2 t 2 



,(17) 



R F 



8tt fl 2 (1 + cos 2 S) 2 + 4w 2 r 2 



(19) 



which is a decreasing function of magnetic field. 

Conclusions - The results of our analysis are summa- 
rized in Fig. 3 which shows the full spin Hall conduc- 
tivity asH, including Rashba, skew-scattering and side- 
jump terms, as a function of two variables, frequency 
and magnetic field, for realistic values of the parameters. 
The figure shows the smooth crossovers between differ- 
ent regimes and should be useful to experimentalists at- 
tempting to extricate the various components of this still 
quite intriguing effect. 
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spin-precession corrections to Eq. are of order of a^ai 
and we omit them as small. 

Substitution of operator for total system Yl =iPiv^i^ by 
sum of one-electron operators Np y S z is justified because 
we omit electron-electron interactions. 
Actually, this is not the complete formula for the skew- 
scattering contribution (which is quite complicated), but 
a simplified version of it, which is valid in the limit 
of OLRkFT ^> 1. However, even in the opposite regime 
oiakFT <C 1 this formula remains qualitatively correct, 
in the sense that it exhibits the correct scaling behaviors 
in the limits of weak and high magnetic field, low and 
high frequency. Notice that contributions of order uj 2 /Q. 2 , 
which we have omitted, give finite SS at finite frequency 
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and zero magnetic field, but are negligible in comparison 
to the larger intrinsic contribution. 
We omitted corrections of order a 2 R cti. 



Again, we present only a simplified version of the complete 
formula, subject to restrictions explained in footnote 18 . 



